Magnetic neutron scattering study of YVO3: Evidence for an orbital Peierls state 
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Neutron spectroscopy has revealed a highly unusual magnetic structure and dynamics in YVO3, 
an insulating pseudocubic perovskite that undergoes a series of temperature induced phase transi- 
tions between states with different spin and orbital ordering patterns. A good description of the 
neutron data is obtained by a theoretical analysis of the spin and orbital correlations of a quasi-one- 
dimensional model. This leads to the tentative identification of one of the phases of YVO3 with the 
"orbital Peierls state" , a theoretically proposed many-body state comprised of orbital singlet bonds. 

PACS numbers: 75.30.Et, 75.30.Ds, 75.50.Ee, 78.70.Nx 



Materials with valence d-electrons exhibit a multitude 
of competing many-body states whose theoretical de- 
scription is still in its infancy. Recently, much atten- 
tion has focused on transition metal oxides with low- 
lying electronic states (termed "orbitals" ) in which tem- 
perature or doping can drive phase transitions involv- 
ing marked redistributions of the valence electron den- 
sity ("orbital ordering") yj. The orbital ordering tem- 
peratures in insulating oxides (such as the widely studied 
manganites) are generally high and approach the temper- 
atures at which these materials become chemically un- 
stable. Magnetic phase transitions involving unpaired 
valence electrons then occur at much lower tempera- 
tures. A striking exception to this general scenario was 
recently discovered in insulating YVO3, where a series 
of temperature-induced magnetization reversals heralds 
a spontaneous redistribution of the valence electron den- 
sity far below the magnetic ordering temperature 0,1110. 

The microscopic origin of this unusual series of phase 
transitions has thus far remained elusive. Given that 
YVO3 is an insulator with a simple, nearly cubic lat- 
tice structure and only two localized valence electrons 
per vanadium atom (3d 2 ), the difficulty of obtaining a 
microscopic description of its phase behavior may seem 
surprising. It is rooted in the large number of nearly de- 
generate many-body states accessible to transition metal 
oxides with unquenched orbital degrees of freedom, a sit- 
uation termed "orbital frustration" . Clearly, the orbital 
occupations of YV0 3 cannot be treated as temperature 
independent parameters and are hence affected by or- 
bital frustration. However, while pure spin models with 
fixed exchange interactions have been extensively stud- 
ied, quantitative solutions of models with variable spins 
and orbitals have thus far been obtained only in one spa- 
tial dimension (ID) In this Letter we show that 



a ID spin-orbital model contains much of the physics 
underlying the phase behavior of YVO3. The tempera- 
ture induced charge density rearrangement is driven by 
a thermal crossover between two competing states near a 
zero-temperature phase transition of the model: a ferro- 
magnet and a novel collective singlet state, the "orbital 
Peierls state", that has been theoretically proposed || 
but hitherto not been observed. In this state, the su- 
perexchange energy of a system of fluctuating exchange 
bonds is lowered through spontaneous formation of a 
dimerized state with alternating strong and weak bonds, 
in analogy to the Peierls instability of a ID metal. 

The neutron diffraction experiments were carried out 
on the E5 four-circle diffractometer at the BER II reactor 
of the Hahn-Meitner Institute in Berlin, Germany, where 
Cu- or PG-monochromators selected the neutron wave 
vectors 7.07 A -1 and 2.66 A -1 , respectively. The data 
were collected with a two-dimensional position sensitive 
3 He-detector. The inelastic neutron scattering experi- 
ments were carried out at the IN22 triple-axis spectrom- 
eter located at the Institut Laue-Langevin in Grenoble, 
France, with a pyrolytic-graphite (PG) monochromator 
and a PG-analyzer, both horizontally focusing. The fi- 
nal wave vector was fixed at either 2.66 A -1 or 4.10 A -1 . 
This allowed the determination of the magnon dispersion 
relations in an energy range up to 50 meV, with energy 
resolutions varying between I and 5 meV. The samples 
were untwinned single crystals of YVO3 (3 x 3 x II mm) 
grown by a floating zone technique described earlier |6j. 

The room temperature crystal structure of YVO3 is 
described by the orthorhombic space group Pbnm 0, 
and lattice parameters a = 5.2821(9) A, b = 5.6144(8) A, 
and c = 7.5283(11) A. For simplicity we will follow com- 
mon practice and index the reciprocal lattice in terms of 
a pseudocubic subcell with (almost equal) lattice param- 
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FIG. 1: Integrated intensities of (a) the (1/2, 1/2, 1/2) and 
(b) the (1/2, 1/2, 0) Bragg reflection of YVO3 as a function 
of temperature. Both reflections are structurally allowed, but 
most of their intensity is of magnetic origin. Pictorial repre- 
sentations of the magnetic structures are given as insets (to- 
gether with the octahedral rotation pattern corresponding to 
the orthorhombic cell). A structural phase transition at 200 
K leads to a weak modification of the structural component 
of the (1/2, 1/2, 0) reflection H0. 



eters a/y/2, b/y/2, and c/2. YVO3 exhibits two magnetic 
phases. Between 1.5 K and Tni = 75 K, we observe mag- 
netic Bragg reflections of the type (h/2, k/2, 1/2) with 
h, k, I odd (Fig. la). The magnetic structure is thus 
of G-type, that is, antiferromagnetic in all three direc- 
tions of the pseudocubic cell (inset in Fig. la). A full 
refinement of the magnetic structure (R(F) = 0.034) 
yields an ordered moment of 1.72(5)/zb oriented along 
c. The magnetic diffraction pattern changes at Tni, and 
for T N1 < T < T N2 = 116 K reflections of the type 
(h/2, k/2, 1), with h, k odd and I = or even, are dom- 
inant (Fig. lb). In this temperature range, the mag- 
netic structure is predominantly of C-type, that is, anti- 
ferromagnetic in the a6-plane and ferromagnetic along c. 
These observations are consistent with previous work on 
powder samples Q. 

However, our single-crystal data now show that the 
high temperature phase is much more complex than pre- 
viously assumed. For instance, magnetic Bragg reflec- 
tions of G-type persist in this phase with diminished in- 
tensity (Fig. la), indicating that the spin structure is 
noncollinear. A refinement of the magnetic diffraction 
pattern (R(F) = 0.093) yields a G-type component of 
0.30(4)^b along c, and C-type components of 0.49(3)^b 
and 0.89(2)/zb along a and b, respectively. The moments 
are thus canted by 9 — 16.5(1.8)° out of the ab plane, 
and the total ordered moment is 1.05(2)/ib, much smaller 
than both the free-ion moment of 2/iB and the ordered 
moment of the low temperature phase at Tni ■ The large 
reduction of the order parameter suggests strong quan- 
tum fluctuations in the high temperature phase. 

The unusual nature of the high temperature phase is 
underscored by its magnetic dynamics. Representative 
inelastic neutron scattering data are shown in Fig. 2, and 
a synopsis of the magnon dispersions extracted from fits 
to these data is given in Fig. 3. The magnons in the low 
temperature phase follow a simple anisotropic Heisenberg 
model with exchange parameters J c — 5.7±0.3 meV and 
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FIG. 2: Solid symbols: Inelastic neutron scattering profiles 
in the high temperature phase of YVO3, showing acoustic and 
optical magnon branches. (Note that I = 1 is equivalent to 
1 = 0.) The lines are the results of fits to a magnon cross 
section convoluted with the spectrometer resolution function 
[T^ |. Open symbols: Data taken in the paramagnetic state. 



J a b = 5.7 ± 0.3 meV (Fig. 3a). The zone-center gaps 
can be attributed to a single-ion anisotropy of the form 
-KiSf z + K[(S? X - Sf y ) with Ki = 0.33 ± 0.08 meV 
and K[ — 0.18 ± 0.05 meV. Up to the first order transi- 
tion at Tni , the G-type magnetic correlations are weakly 
temperature dependent (Fig. la). In the high temper- 
ature phase, low energy magnons are observed at both 
Q = (1/2,1/2,1/2) and (1/2,1/2,0), as expected based 
on the canted spin structure (Fig. lb). The magnon gaps 
at these points are identical and substantially smaller 
than the zone-center gap of the G-phase immediately be- 
low Tni, ruling out any two-phase coexistence scenario. 

The spectrum of the high temperature phase (Figs. 
2 and 3b) exhibits several unexpected features. First, 
the magnon band width along the ferromagnetic c-axis 
is larger than that in the antiferromagnetic afe-plane. 
This violates the standard Goodenough-Kanamori rules 
according to which ferromagnetic superexchange inter- 
actions are generally substantially weaker than antifer- 
romagnetic interactions. Further, the spectrum is split 
into optical and acoustic magnons with a gap of 5 meV 
between both branches. The optical/acoustic splitting, 
which constitutes a large fraction of the total magnon 
band width, can only be described by assuming two dif- 
ferent ferromagnetic exchange bonds along the c-axis, so 
that the spin Hamiltonian becomes 



h= J *»& ■ 3o - E K ^ - E w ± s ) 

(i,j)eab i (ij')ec 

{Si ■ Sj) KzSixSjx i d(Si X Sj Z Si Z Sj X )] (1) 

where S x \\ c, K\ and K2 are the single- ion and symmet- 
ric superexchange anisotropies, and the Dzyaloshinskii- 
Moriya vector d is along the octahedral tilt axis which is 
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FIG. 3: Magnon dispersion relations in the (a) low and (b) 
high temperature magnetic phases. The lines result from di- 
agonalizations of anisotropic Heisenberg Hamiltonians in lin- 
ear spin wave theory, as discussed in the text. The effective 
line widths, indicated by the bars, can be attributed to the 
superposition of two unresolved magnon branches. 



staggered along c. An excellent fit of the spectrum can 
be obtained by J ab = 2.6 ±0.2 meV, J c = 3.1 ±0.2 meV, 
S = 0.35, K x = 0.90 ± 0.1 meV, K 2 = 0.97 ± 0.1 meV, 
and d — 1.15 ± 0.1 meV (Fig. 3b). The canting angle 
extracted from Eq. 1,6= \ tan" 1 [2d/(2J c -K 1 -K 2 )] 
1 - 1 ■ 1°, also agrees well with the observed value. 

Because of the crystallographically indistinguishable 
V-V distances and bond angles along the c-axis p|, our 
observation of strongly alternating exchange bonds in 
this direction must be regarded as a signature of highly 
unusual orbital correlations. We now proceed to outline 
a quantum many-body description of these correlations 
and their influence on the spin dynamics; details will be 
published elsewhere. According to electronic structure 
calculations 0, , the vanadium ^g-levels are split into 
a lower-lying singlet of xy-symmetry and a higher-lying 
doublet spanned by the xz- and yz-orbitals. A structural 
transition associated with a contraction of the c-axis Q 
indicates that this basic hierarchy of electronic states 
is established around 200 K. The strong intra-atomic 
Hund's rule interaction stabilizes a high-spin (S = 1) 
state in which the half-occupied xy-orbital dominates 
the antiferromagnetic in-plane interaction at all tempera- 
tures below 200 K. Below the structural phase transition 
at 77 K, the degeneracy of all three i2 9 -levels is split by 
lattice distortions, resulting in rigid orbital order with 
antiferromagnetic exchange coupling in all directions as 
observed experimentally. 

According to a recent theory |l0j ]. however, an unusual 
electronic state with orbital order involving only the xy- 
orbitals is established above 77 K. In this state, the de- 
generacy of the xz- and yz-orbitals, which control the 
superexchange interaction along the c-axis, is not lifted 
by lattice distortions. Coupling to spin excitations via 
the superexchange Hamiltonian broadens the manifold 
of degenerate quantum states spanned by these orbitals 
into a band of overall width J = At 2 /U (where t is the 
V-V hopping parameter and U the intra-atomic Coulomb 
interaction). The band of correlated spin-orbital fluctu- 
ations is one- dimensional, because only c-axis bonds are 
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FIG. 4: (a) Nearest-neighbor pair spin correlation function 
{Si ■ Si+i) of the ID spin-orbital model described in the text 
for A = and different values of w — Jh/U, and of a ID fer- 
romagnet with fixed exchange interactions, as a function of 
temperature, (b) Correlation length £ (corrected for the lead- 
ing 1/T divergence for critical excitations) for 77 = 0.15, and 
A = and A = 0.1J respectively. For the latter parameter 
set the spin anisotropy (S^) — | is also shown. 



involved. (Spin and orbital interchain interactions are 
more than an order of magnitude smaller than J.) We 
now show that this band is unstable against dimeriza- 
tion, and that this "orbital Peierls" instability holds the 
key to the unusual magnetic properties we have observed 
in the C-type phase. 

In addition to the overall superexchange energy scale 
J, the spin-orbital Hamiltonian along the c-axis contains 
two parameters: the Hund's rule coupling Jh (via the 
ratio 77 = Jh/U), and the spin orbit (SO) coupling A. 
To develop some intuition, we first consider the influence 
of the largest parameter, J, separately by setting n = 
A = 0. In this limit, the Hamiltonian can be written |l(| 



H =o £(£-^ + l)(7 



+ 7) 



(2) 



where < ij > denotes a nearest- neighbor pair of V-atoms 
along c, and f is a pseudospin-1/2 operator acting in the 
xz- and yz-orbital subspace. The energy of a single ex- 
change bond described by this Hamiltonian is minimized 
by formation of an orbital singlet ((t^ •Ti+i) = —3/4) with 
an accumulation of charge between the V-atoms and fer- 
romagnetic alignment of the spins. Hence a strong dimer- 
ization in the ground state is expected [ll| . 

In order to elucidate the ground state properties of 
the full many-body Hamiltonian, we performed numeri- 
cal calculations using the transfer matrix renormalization 
group (TMRG) method which is known to yield accurate 
finite-temperature correlation functions and thermody- 
namic properties for a variety of ID systems 
Some results are shown in Fig. 4. The numerically com- 
puted correlator (fi • f^+i) extrapolates to -3/8 per bond 
at T=0, thus demonstrating that the ground state in- 
deed spontaneously dimerizes into orbital singlets on ev- 
ery second bond as suggested by the simple consideration 
above (data not shown). While the exchange interaction 
between spins within an orbital singlet is strong and fer- 
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romagnetic, nearest-neighbor spins in adjacent singlets 
are coupled through a weak, antiferromagnetic exchange 
interaction; the spin pair correlator (Si -S i+ i) is therefore 
negative (Fig. 4a). 

The spontaneous exchange bond dimcrization in the 
ground state of the Hamiltonian Eq. [21 already hints 
at an explanation of the unusual magnetic spectrum ob- 
served in the high temperature phase of YVO3. How- 
ever, a quantitative description of YVO3 must also in- 
clude the Hund's rule and SO interactions. The gener- 
alization of Eq. El for nonzero Hund's rule coupling is 
given in Refs. |lOt fl3j|. and the numerically computed 
correlator (Si ■ Si+i) is plotted in Fig. 4a for different 
values of r\. As expected on general grounds, the ground 
state for large r] becomes a uniform ferromagnet with 
(Si-Si+i) = 1. The transition between dimerized and fer- 
romagnetic states occurs at rh. ~ 0.11, remarkably close 
to the realistic value ~ 0.12 |9| for vanadium oxides. 

Due to the proximity to the zero-temperature phase 
transition, the spin and orbital correlations at elevated 
temperatures exhibit an intriguing evolution. First, a 
pronounced maximum of the spin correlator (Si ■ 
as a function of temperature is apparent for the case of 
77 < 77c . The origin of the initial dramatic increase upon 
heating lies in the weak exchange bonds between spins in 
adjacent orbital dimers. These bonds are highly suscep- 
tible to thermal fluctuations, and the associated increase 
in entropy stabilizes the orbital Peierls state. In fact, 
upon increasing the temperature this state quickly be- 
comes the leading instability even for 77 > r/ c , and the 
numerically computed dimer correlation length begins to 
exceed those of other possible ground states (Fig. 4b). 
With J ~ J a b/0.16 ~ 40 meV extracted from the low 
temperature magnon dispersions following Ref. 01 > the 
temperature range in which the dimerized phase is the 
leading instability coincides approximately with the high 
temperature phase of YVO3. This indicates a micro- 
scopic mechanism for the exchange bond dimerization ex- 
perimentally observed in this phase. Fig. 4a also shows 
that the spin correlations of the spin-orbital chain are 
frustrated and hence weaker than those of an orbitally 
nondegenerate ferromagnetic chain. The suppression of 
the spin correlations due to orbital fluctuations provides 
a natural explanation for the anomalously small ordered 
moment observed in the high temperature phase. 

Next, we address the intra-atomic SO coupling Hso = 
—2ASi Z Ti Z where S z || c. TMRG calculations in which 
this interaction is explicitly included (Fig. 4b) show 
that the interplay between spin-orbit and superexchange 
interactions gives rise to an easy-plane spin anisotropy 
while the entropy-driven dimcrization is almost unaf- 
fected. In the temperature range of interest, the quantity 
(S%) drops below |, its value for an isotropic paramagnet. 
The preferred spin direction is thus the afr-plane, as ex- 
perimentally observed in the high temperature phase of 
YVO3. These nonlocal orbital correlations are at the root 



of the large superexchange anisotropies of Eq. 1. They 
counteract the single-ion anisotropy (which is responsible 
for the c-axis oriented moments in the low temperature 
phase) and account for the noncollinear spin structure 
in the high temperature phase. The physical origin of 
this effect is a compromise between the tendency of the 
superexchange interaction to co-align the spins and anti- 
align the orbital pseudospins of neighboring V-atoms, 
and the tendency of the SO interaction to align spin and 
pscudospin at every site. In analogy to a spin- flop transi- 
tion (with the c-axis oriented orbital pseudospin playing 
the role of the magnetic field), a ferromagnetic alignment 
of a6-oriented spins satisfies the superexchange coupling 
while allowing the spins to cant out of the plane to take 
partial advantage of the SO coupling. 

In summary, all of the unusual features of the high 
temperature phase of YVO3 including the large ferro- 
magnetic exchange coupling along the c-axis, the optical- 
acoustic splitting of the magnon spectrum, as well as the 
easy-plane anisotropy, the anomalously small magnitude 
and the large canting angle of the ordered moment are 
thus explained in a quasi- ID model with dynamical spin 
and orbital degrees of freedom. Based on this analysis a 
lattice dimerization is expected in the high temperature 
phase, but may be difficult to observe because of the large 
thermal and quantum fluctuations in the orbital sector 
and the weak lattice coupling of the ti g orbitals. 
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